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Abstract 

We apply the renormalization group approach to the sandpile on the triangular 
lattice. The only attractive fixed point is found. The obtained fixed point height 
probabilities are compared with numerical simulations. The value of critical expo- 
nent of avalanche size distribution is found to be r = 1.36. The probabilities of the 
sand transition are compared with the branching probabilities of the spanning trees 
on the triangular lattice which are also evaluated. 
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1 Introduction. 



The main idea of the renormahzation group (RG) approach to critical fenomena is the 
calculation of parameters of a critical state by a successive transition from small-scale 
interactions to the large-scale ones. It is possible due to the vanishing of characteristic 
temporal and spatial scales at the critical point leading to temporal and spatial self- 
similarity. A wide range of the systems exhibiting these properties is unified by the 
concept of Self-Organized Criticality (SOC) introduced by Bak et a/ These systems 
are characterized by the evolution into the critical state free from fine tuning of any 
parameters. The sandpile model was proposed in the same work to illustrate the idea of 
the SOC. This cellular automaton was found to be rather a simple object of investigations, 
although rich in content. The sandpile model involves the set of sites on which the function 
of height is defined. Time after time being added to a random site, "grains of sand" cause 
the growth of its height. If the height reaches some critical value, the site topples reducing 
its height and increasing the heights of several neighboring sites. In turn, if the heights 
of the neighboring sites reach the critical value, they also topple and so on. The "sand" 
can leave the system on the boundary. After evolving for some time the system comes 
into the recurrent subset of all its configurations which is the critical state. Further the 
successive topplings between the stable states called avalanche return the system into the 
critical state each time it leaves it. 

The most successful RG approach to the description of the sandpile on the square 
lattice has been proposed by Pietronero a/ 1|] . It is based on the real space RG scheme 
driven by steady state condition feedback. The real space renormahzation is performed 
according to the standard block-transformation scheme. The feedback mechanism allows 
one to take into account the strong correlations appearing at the critical point. Practically, 
this approach identifies the dynamic and static parameters of the model. While the 
static configurations are described in terms of stable and critical cell distribution, the 
dynamic renormalized variables characterize the toppling rule for the critical cells. The 
conservation law valid in the critical state makes it possible to obtain the feedback relation 
linking critical values of the dynamic parameters and the static configuration weights. 
The obtained renormahzation equation has a unique attractive nontrivial fixed point. By 
considering the avalanche as a relaxation process at the definite scale the exponent r of 
the distribution of the avalanche size was expressed through the fixed point parameters. 

The next step in this direction was made by Ivashkevich [Q. While previous approach 
dealt only with the critical cell distribution, the cells with different heights were introduced 
into this model. Further, the sandpile evolution was identified with the branching process 
and the kinetic equations for the sand transfer were written as equations of chemical 
reactions of a special type. This finer consideration allows one to take into account more 
processes of toppling inside a cell. In addition, the generating function method proposed 
in the same work makes all calculations much simpler. Thus, the exponent r was defined 
with higher accuracy, and different height site concentrations were calculated. 

Due to the fixed point independence from the smallest scale dynamics, all the results 
are valid for the model with an arbitrary rule of toppling. Therefore, this is correct for 
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the Abelian sandpile model (ASM), which was shown to be exactly soluble by Dhar 
In fact, the values of fixed point concentrations are, indeed, in good agreement with the 
exact values calculated by Priezzhev 0, and value r obtained in and is very close 
to T = 1.25 proposed by Priezzhev et a/ in [0. 

The relationship between the configurations of the ASM and spanning trees was found 
in [^] for the first time. Later the one-to-one correspondence between the allowed ASM 
configurations and the set of spanning trees was shown in 0. In the letter 0],one inter- 
esting fact has also been noted: the fixed point probabilities of the sand being transferred 
from the toppling cell to one, two, three or four neighboring cells were found to be al- 
most equal to the "branching probabilities of spanning trees" which are the probabilities 
of a random site of the spanning tree having a coordination number 1,2,3 or 4, exactly 
calculated in fl^.It makes sense to verify the hypothesis about the equivalence of these 
probabilities and in this way to clarify a possible relation between the RG approach and 
the spanning tree statistics. 

Excellent results of the RG approach on the square lattice call for obtaining the same 
results on other lattices. In our paper, we consider the RG approach for the triangular 
lattice. We compare the obtained equilibrium concentrations with the computer simu- 
lations and the exponent r of the avalanche size distribution with its exact value. We 
also calculate the branching probabilities of the spanning trees on the triangular lattice 
and compare them with the fixed point probabilities of the sand transfer into neighboring 
sites. It is shown that the hypothesis about their coincidence at least does not have a 
universal character. 



2 Renormalization approach. 

Let us consider a sandpile on the triangular lattice. The height of each site can take an 
integer value Z = 1,2, ... . If the height of any site exceeds the critical value Zc = Q, it 
topples. At this moment, the particles can be transferred to one or several neighboring 
sites, increasing their heights by one. This toppling rule is characterized by the vector 
P = {PiiP2, ■■■yPe), where pt is the probability for the site to topple into i neighboring 
cells. 

Once we want to bring the consideration of this model from the smallest scale to a 
large one, we should define the successive scale increasing procedure. The step of this 
procedure consists in the transition from the sublattice Lb, which is constructed of the 
cells of size b on the initial lattice L, to the sublattice -Z^^Vs "^ith the cells of size 6^3. To 
realize it, three cells of the sublattice Lb are unified into the block. This block is replaced 
by the cell of Lj^^, as is shown in Fig. la. All bonds of Lb outgoing from the given block 
and coming into another block are included into one renormalized bond of L^^^, as is 
shown in Fig lb. 

Now, we must describe the coarse grained sandpile dynamics. In the following calcu- 
lations we use the RG approach in the form proposed in and 0]. Let us define the 
coarse grained variables for an arbitrary scale. They are introduced in such a way that 
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one can treat the sandpile properties in any scale like ones of the original sandpile. Each 
cell of Lb is characterized by two vectors. The first one is n*^'') = (n^, ub, nc, ud, ue, np)- 
This is the set of probabihties for the Lf, cell to behave, with respect to transformations 
initiated by the dropped into the cell particle, like a site on the initial lattice with height 
1,2,3,4,5 or 6, respectively. In particular, np is the probability that the cell is critical. 
The second vector, p*^*) is the analogue of the vector p in the initial sandpile. 

Let us represent the coarse grained dynamics of the sandpile as the following branching 
process. 
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Here the species A, B, C, D, E, F are the cells which behave like the sites with heights 
1,2,3,4,5,6 , and (p,(p are the particles obtained by the cell and transferred into other 
cells, respectively. This branching process can be reinterpreted as irreversible chemical 
reactions for which we can write the kinetic equations of the transport of " sand'' : 



flA = 






flB = 




{P5 np + nA - ub) 


nc = 




{P4 uf + nB- nc) 






(P3 np + nc- nj)) 


riE = 




{P2 np + riD - ns) 






(pi np + ns - np) 






(p np - 1) + ^ z/V^(n^np) + ?7(r, t) 


p = 


6pl 


+ hpl + Apl + 2>pl + 2pl+p\ 



(2) 



Here n^ denotes the concentrations of particles ?](r, t) is the noise term corresponding 
to the randomly dropping particles, f is the diffusion coefficient ,which is equal to ^ for 
the triangle lattice, p is the average number of particles outgoing from the toppling cell. 

Performing the renormalization procedure we express the parameters of through 
ones of Lfo. 
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The critical values of p** and obtained by taking the limit b oo satisfy the renormal- 
ization equations 

n* = n(n*,p*), p* = p(n*,p*) (4) 

Due to nonlocal properties of the sandpile dynamics, writing down the first renormaliza- 
tion equation does not seem possible. To avoid it, we can exclude n* from the second one. 
We use equations (2) and the requirement for the SOC state to be stationary. Taking 



n* = 0, (5) 



we can directly obtain 
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n*^ = P*^+P*^ + P*^+P^ (6) 
P 
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Pl+P*5+Pl+P3+ Pi 
P 

pI + pI + pI + pI + p*2+ pI ^ 1 
P p 



Substituting this into the second equation we obtain a closed renormalization equation 
for the vector p. Search for the relationship between p^^^ and 

p{bV3) is 

reduced to counting 

the processes at the scale b leading to the toppling of the cell at the scale fey^. To this 
end, we, following introduce the generating function of the toppling of the cell at 
an arbitrary scale. If we denote six directions outgoing from the cell by xi, ...,X6, the 
generating function of toppling of the sublattice cell will be given by the formula 

P? 

a{xi,X2, X3, X4, X^, Xe) = ^{Xi + X2 + X^ + X^ + X^ + Xq) + 



-Y^{xiX2 + X1X3 + X1X4 + ... + x^xe) + -^(^" + (7) 

— — cr"" + ——a'"" + pe^''>{xiX2XsX4X5Xe), 
15 o 

where a with i primes is the sum of all possible products Xk^Xk2---Xki of i different terms 
Xfc^, ■■■iXki ■ This function counts all processes leading to sand transition from the cell of 
Lfc. The coefficient of each polynomial term is the probability for particles to go to the 
corresponding directions after the cell toppling. This function has the following properties: 
a) If the arguments corresponding to any directions are replaced by zero, the function 
counts the toppling processes that don't send the particles to these directions. 
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b) If they are replaced by unit, the function counts the topphng processes irrespective 
of outgoing to these directions. 

c) If we want to describe the toppling of a couple of cells, where the particle sent 
by the first cell induces the toppling of the second one, we must replace the argument 
labehng the direction of outgoing of this particle by the a-function of toppling of the 
second cell. This function contains also the processes in which the particle does not go in 
the given direction and the second cell does not topple. It is necessary to subtract them, 
if we want to take into account only processes spreading to both cells. 

d) In addition, it is easy to see that the function is normalized so that cr(l, 1) = 1. 
Combining a-functions in different ways we can build the functions that count the 

specific topplings. Once we manage to link the generating functions at different scales, we 
express the probabilities p at a given scale through ones at a new scale. To this end, we 
write the generating function for the block consisting of three cells of Lf,. This function 
must count the processes leading to the transition of particles from this block to some 
neighboring blocks. We can perform this using cr-function properties outlined above. In 
addition, according to [0] we must consider only processes which match the spanning 
condition, i.e. with the characteristic size equal to the size of a cell. In our case, these 
are relaxation processes in which more than one site of Lf, participates. All three types 
of blocks, whose topplings match the spanning condition and their relaxation schemes 
are shown in Fig 2. The rest of the schemes can be obtained from the shown ones by 
rotations. 

For going over to the sublattice L^^, we replace the block considered by the cell of 
Lfj^, respectively renaming the outgoing directions, as is shown in Fig. lb. To give the 
obtained function the meaning of the generating function, we must subtract all processes 
not going out of the block and normalize the result. Thus, we obtain the generating 
functions y2, ys, y^, j/s, y&){3 = a, b, c) for the topplings of blocks of three kinds: a, b 

and c. For example, the generating function for relaxation of the block of type a is the 
following: 



Sa(yi,2/2,?/3,y4,y5,2/6) 



^ajyi, ■■■,1/6) - Sa(0, ...,0) 

S,(1,...,1)-S,(0,...,0) 



where 



Sa(?/l,?/2,Z/3, 2/4, 2/5,1/6) = 

(0-(2/l, 0^(^2,2/2, 2/3,1,1,2/1), 1,2/5, 2/6, 2/6) - O-(yi,0, 1,7/5,2/6,2/6) + (9) 

(7(2/2, 2/2, 2/3, l,cr(2/i, 1,1, 2/5, 2/6, 2/6), 2/1) - o-(2/2, 2/2, 2/3, 1,0,2/1)) + c.p. 

The arguments yi ...ye denote the directions on the lattice at the new scale. The terms c.p. 
are obtained from the first term by the cyclic permutations of arguments corresponding 
to the rotation by the angles ^ and The generating functions for the blocks of 
other types look analogously but they are too cumbersome to be given here. To write a 
general generating function, we accomplish summation over all types of blocks. Different 
types of blocks appear with the different probabilities which depend only on their heights. 
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Therefore, we must multiply their generating functions by the corresponding weights. 
After summation the result must again be normalized 



2/2,2/3, 2/4, ?/5, ye) 



^(yi,-,y6) 

E(l,...,l) 



(10) 



Wa = 3nF^(nA + nB + nc + no), = Suf^ue, Wc = np^ 

Thus, we have the normalized generating function describing the toppling of a cell of the 
lattice -Z^bVs- According to the RG ideology, all particles outgoing from a cell in some 
direction become one renormalized particle. Therefore, the second powers of yi must be 
changed to the first powers. Now the coefficient of each term is again the probability of 
the respective process of toppling of the size cell. Taking the coefficients of sums 
analogous to (1), we obtain new probabilities p(*^) expressed through p^^^ and n^''\ 
Eventually, we have sought a system of renormalization equations. Solving this system, 
we obtain the fixed point values of p and n.They are shown in tables 1 and 2. 

According to Pietronero et al, the critical exponent of the avalanche distribution r 
may be expressed through the critical point parameters as follows. The probability of 
an avalanche of the linear size r is P(r) dr = r*^^"^"^^ dr. Thus, the probability for the 
avalanche to cover the cell of the size b dying beyond it is given by the formula: 

^ - TO? - ' - (^)^"-^' 

On the other hand, K can be expressed through the fixed point parameters 

K = pI{1 - np) + plil - n*pf + ... + pl{l - n^^f (12) 
Using these expressions we obtain 

r = 1 - MIH^ ^ 1.366 (13) 
In 3 ^ ^ 



3 Results and discussion. 

Obtaining the height probabilities in the sandpile on the triangular lattice by the RG 
approach seems to be efficient because of the absence of an exact result. While the RG 
approach is relatively simple, the exact calculations require applying a very complex tech- 
nique even in the simplest case. In Table 1, we compare the fixed point height probabilities 
with the results of numerical simulation .One can notice a very similar behavior of the 
numerical and RG results. 
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n*2 


nl 


nl 


nl 


nl 


RG 

Numerical 


0.036 
0.054 


0.135 
0.092 


0.197 
0.139 


0.209 
0.188 


0.210 
0.240 


0.210 
0.281 



Table 1. - Fixed point height probabihties in comparison with the numerical results. 
The numerical simulations were performed on the lattice 300 x 300 with statistics 10^ 

configurations. 

Some distinction of their values is stipulated by initial roughness of the RG approach 
concerned with discrete increasing of the scale and averaged consideration of the coarse 
grained dynamics. However, the advantage of this RG scheme is a quite attractive char- 
acter of the obtained fixed point. 

The obtained value r = 1.366 is close to the exact one 1.25 although the coincidence 
is rather worse than in the square lattice. This can be explained by a small size of a 
renormalized cell. The renormalization scheme thus constructed takes into account only 
a minimal number of toppling processes. Therefore, these calculations don't give high 
accuracy. However, this result allows us to speak about the stability for a given RG. 

At the end, in Table 2, we compare the fixed point probabilities of the sand transfer 
Pk with the branching probabilities of spanning trees calculated in Appendix. 
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0.0000179 


0.00226 


0.0558 


0.296 


0.471 


0.174 


Branching 














probabilities 


0.322 


0.417 


0.207 


0.0488 


0.00553 


0.000241 



Table 2. — Fixed point probabilities pk of the sand being transported to k neighboring cells 
and branching probabilities of spanning trees on the triangular lattice. 

The comparison of the probabilities of the sand transfer with the branching probabil- 
ities of spanning trees is very surprizing. As these quantities for the square lattice are 
very close, a similar behavior can be expected in our case as well. However, one can see 
quite an inverse order of the obtained values for the example considered. 

Thus, we can conclude that the hypothesis about their coincidence is not confirmed 
in the case of triangular lattice. Further clarification of this question seems to be related 
to more complete investigations of a sandpile on a square lattice. 
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Appendix 



The branching probabihty of a spanning tree pk is the probabihty for any site of a random 
spanning tree to have a coordination number k. The effective tool used for evaluating 
branching probabilities is provided by the Kirchhoff theorem about spanning trees. The 
method of applying it is not novel. Below we follow the methods outlined in the article 
1^. Here we recall its main features. To formulate the Kirchhoff theorem we consider the 
arbitrarily connected graph L consisting of n + 1 sites. Let us fix any site which will be a 
root. The weight Xij corresponds to each bond of L edged by the sites i and j. Let A{x) 
he a n X n matrix: 

Aij{x) = ^ —Xij , i and j are adjacent sites (14) 
, otherwise. 

While the summation is performed over all sites of G i and j run over all sites different 
from the root. Then, according to the Kirchhoff theorem, 

g{x) = detA{x) (15) 

is nothing but a generating function of the rooted spanning trees on this graph. If all Xij 
are equal to 1, the matrix 

A = A{x,, = 1) (16) 
is the discrete Laplacian on graph L and 

N = detA (17) 

gives a total number of the spanning trees on this graph. 

For calculating a number of spanning trees on the triangular lattice L whose random 
site i has k adjacent sites, let us consider the lattices L' with the defect in the cell 
containing site i. In other words, we deal with lattices obtained from L by deleting one 
or several bonds outgoing from i. If we leave only one bond connected with i, all trees 
on this new lattice will occupy it. Summing over all possible directions of this bond, we 
obtain the number of spanning trees the site i of which has only one adjacent site. In the 
same manner we can link the numbers of spanning trees with other coordination numbers 
of site i with total numbers of spanning trees on the lattices with a local defect 

= E (18) 
fe=i "-^e 

where A''^ is the number of spanning trees in which the site i has the coordination 
number k. N'j^ is the sum of total numbers of spanning trees on the lattices obtained from 
L by deleting all bonds connected with site i except k ones which is taken over all possible 
positions of these bonds. Dividing A^^ by the total number of spanning trees, we obtain 
the branching probabilities. 
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Thus, we reduce the problem of searching for branching probabihties search to the 
calculation of the numbers of spanning trees on the lattices with the defect. To this end, 
we use the Kirchhoff theorem. According to it, the number of spanning trees on the lattice 
with a defect is given by the formula 

N' = detA' (19) 

A' is the discrete Laplacian on a new lattice. If we refer to the difference between A and 
A' as a defect matrix 6, branching probabilities are expressed through the terms of the 
following form: 

TV' detA' 

/ is the unit matrix; G satisfies the matrix equation 

AG^I (21) 

In the thermodynamic limit G is nothing but the Green function of the Laplace equation 
on the triangle lattice. For an infinite lattice G depends only on the difference of site 
coordinates rather than their values 

q(y - r ) = -^2) +/?(?/i -1/2)) - 1 da d^ 

Jo Jo 1 - l{cos a + cos/3 + cos{a + p)) 27r 2-K ^ ' 

Due to defect localization, the matrix 5 has only a finite number of nonzero elements. 
Thus, we should calculate the 7x7 matrix determinants. In this way, we obtain the 
following results: 

25 324^3 540 66^3 14 55 

Pi = + —-R -r + —5- + - 



108 7r5 TT^ 12v^7r 

235 1620^3 3213 702^3 277 223 
P2 = -TTT -7 — + + 77-^ + 



144 7r5 TT^ 7r3 2 7r2 12 



17 3240^3 7452 2160^3 838 241 
P3 = ^ + -5 - + 3 ^ + TTir- 23 

36 TT^ TT"^ TT^ 2 V 3 TT 

2489 3240 \/3 8478 2928^3 1489 2203 



216 7r5 ^ ^ tt^ 6^7r 

115 1620 \/3 4752 1854^3 1088 3857 
P5^ — 7r + — -R + — -5 + 



9 TT^ n'^ 12y3 7r 

175 324^3 1053 450^3 569 355 
Pe^-TTT -F— + —7 + 



48 7r5 TT^ TT^ 2 7r2 4\/37r 
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Figure Captions 



Fig.l: a) Transformation from the lattice Lh to L^^^. It is performed in such a way 
that the new lattice is again triangular. 

b) Block of three cells of the sublattice L^. This block becomes a cell on the 
new lattice. We show that the directions outgoing from the cells, forming the block, 

(2) (2) 

are coupled to the directions of the lattice at the next scale. For example: x\ ,X2 

(2) (3) 

are replaced by yi and ,X2 are replaced by y2- The superscript denotes the cell 
from which the given direction outgoes. 

Fig. 2: Three types of blocks, whose toppling should be taken into account, and their 
relaxation schemes. The other blocks can be obtained from these ones by rotations. 
The letters in the circles denote the height of the cells. X can take the value A,B,C 
or D. 
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(b) 
Fig.l. 
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(b) 



(c) 




Fig. 2. 
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